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The lingering problem with high central densities in dark halos has arisen in the context of 
(L)CDM cosmologies with n = 1 scale-invariant initial power spectra. Although n = 1 is often 
justified by appealing to the inflation scenario, the choice is not generally justified. Specifically, 
inflation models with mild but important deviations from scale invariance (n ~ 0.9) are not uncom- 
mon, and those with significant "running" of the spectral index are quite plausible. Even a mild 
deviation from scale invariance can be important because halo collapse times and densities depend 
on the relative amount of small-scale power. Here, we choose several popular, often well-motivated, 
models of infiation and work out the ramifications for galaxy central densities. For each model, we 
calculate its COBE-normalized primordial power spectrum and deduce the implied halo densities 
using a semi-analytic method calibrated against N-body simulations. We compare our predictions 
to a sample of ~ 50 dark matter-dominated galaxies using a non-parametric measure of the density, 
Av/2, defined as the mean mass density, relative to the critical density, within the radius at which 
the rotation curve falls to half of its maximum value. While standard n = 1 LCDM halos are 
overdense by a factor of ~ 6, several of our example inflation-|-CDM models predict halo densities 
well within, and even below, the range preferred by observations. We also show how the presence 
of massive (m^ 0.5 eV) neutrinos can help to alleviate the central density problem, even with a 
scale invariant spectrum. We conclude that galaxy central densities may not be as problematic for 
the CDM paradigm as is sometimes assumed: rather than telling us something about the nature of 
dark matter, galaxy rotation curves may be telling us something about infiation and/or neutrinos. 
An important test of this idea will be an eventual consensus on the value of as, the rms overdensity 
on the scale 8/i~^Mpc. Our successful models tend to have values of as ~ 0.75, which is well within 
the range of recent determinations. Finally, models with n > 1 (or as^ 1) are highly disfavored. 

PACS numbers: 98.80.-k, 98.35.Gi, 98.62.Gq, 98.80.Es, 98.80.Cq, 14.60.Pq 



I. INTRODUCTION 

The standard model of structure formation (LCDM) 
is one in which the universe is dominated by cold, col- 
lisionless dark matter (CDM), made flat by a cosmo- 
logical constant (A), and endowed with initial density 
perturbations via quantum fluctuations during inflation 
101 . Although the need for the cosmological constant 
component is unexpected, the CDM-|-inflation paradigm 
is strongly motivated, and with the parameter choices 
f^M = 1 - f^A « 0.3 - 0.5, h « 0.7, LCDM can ac- 
count for an impressive range of astronomical observa- 
tions on large scales. However, on galactic scales, this 
model faces some potentially ruinous difficulties. Perhaps 
the most troubling of these problems are the indications 
that the central regions of galaxies are noticeably less 
dense than the favored LCDM model typically predicts 
|, |, |, |, |, 0, I |, 0. In this paper we suggest that 
these observations might be telling us something funda- 
mental about the epoch of inflation, and also explore how 
massive neutrinos affect halo densities. 

There are two related but distinct facets of the central 
density problem. The first problem concerns the fact that 
the integrated mass densities within well-defined central 
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radn of observed galaxies (see Sec. IV B) seem to be a 
factor of ~ 6 — 8 larger than the central densities pre- 
dicted by standard LCDM H § |. The second, often 
referred to as the "cuspy halo problem," highlights the 
fact that CDM halo density profiles are predicted to di- 
verge at small r (p cx r~",a ~ 0.7 — 1.5), while galaxy 
rotation curves are often better fit with constant density 
cores |, |, |, ||, 0. While the second of these issues 
has received the most attention, it is often degenerate 
with the first problem and is somewhat more controver- 
sial. For example, de Blok and collaborators (lO| have 
argued that their sample of low-surface brightness galax- 
ies favor fits to density profiles with a constant density 
central "core" over those with cusps; however, van den 
Bosch and Swaters argue that a majority of galactic ro- 
tation curves are acceptably fit by divergent profiles as 
long as they are much less centrally concentrated than 
typical halos in standard LCDM ||ll[. Furthermore, all 
observational errors {e.g., slit offset) tend to favor con- 
stant apparent central density over cusps. At present, it 
is not clear that the cuspy halo problem presents a se- 
rious challenge to LCDM, although it appears that the 
data do prefer halos that are less centrally concentrated 
than typical halos in the standard LCDM model. 

The problems with central densities have triggered a 
growing concern that we are missing something funda- 
mental in our understanding of galaxy formation. (This 
is in spite of the fact that some of the problematic 
claims are disputed ||ll|, ||, |l|, |l|l.) Solutions to 
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the problems range from those that use baryonic physics 
111 to those that rely on alter- 

ing the nature of the dark matter ||l|, |2|, |2j, p5| . 
While the astrophysical solutions are reasonably well- 
motivated, the fact that problems seem to exist for rel- 
atively small, dark matter-dominated dwarf galaxies all 
the way to large ellipticals suggests that a single baryonic 
solution may not be able to address all of our concerns. 
The altered dark matter solutions, on the other hand, 
could possibly be made to match the range of observa- 
tions, but only by invoking unmotivated or fine-tuned 
candidates (there is no well-motivated warm dark mat- 
ter candidate). 

Our inflation-derived solution is nonastrophysical but 
works entirely within the desirable tenets of the CDM 
paradigm. This paper is conservative in that it concen- 
trates on the integrated density within some well-defined 
radius, which is certainly more robustly determined in 
simulated halos than the central slope of the density pro- 
file. As we shall discuss below, the same is likely true for 
observed galaxies. 

Our work is principally inspired by Alam, Bullock, and 
Weinberg |^ who suggested that the central density prob- 
lem would be reduced significantly in LCDM if the ini- 
tial inflationary power spectrum were tilted to favor large 
scales. The term "tilted" is defined in terms of the pri- 
mordial power spectrum of density fluctuations, which 
we assume over some range in wave number k can be 
written as P{k) oc fc", corresponding to a mass variance 
per logarithmic interval of A^(fc) = k^P{k)/2n'^. Tilted 
power spectra refer to those with n ^ I. In "standard" 
LCDM, it is assumed that n is exactly 1, corresponding to 
a scale invariant, Harrison-Zel'dovich primordial power 
spectrum. This choice is often justified by the tendency 
for inflation models to predict nearly scale invariant spec- 
tra; however, generic models of inflation do not predict 
primordial power spectra that are exactly scale invariant. 
As we discuss below, the central densities of dark ha- 
los are extremely sensitive to the amount of small-scale 
power and hence small deviations from scale invariance 
can be very important. 

Our aim is to compute the predicted primordial P{k) 
for several example inflation models and apply these re- 
sults to the question of galaxy central densities and con- 
centrations. Although similar in spirit to the agenda of 
Kamionkowski and Liddle who suggested that the 
small-scale crisis might reflect a s harp feature in the infla- 
tionary power spectrum (see Sec. H C), our mind-set is to 
look at models that are not particularly fine-tuned. We 
simply choose fairly representative, simple, single field 
inflationary models and we examine a range of predicted 
power spectra, even one with with a "blue" initial spec- 
trum (n > I). In the context of slow roll inflation, mod- 
els that predict significant tilt generally yield effective 
spectral indices that are strongly scale-dependent or ex- 
hibit significant "running" of the spectral index. Con- 
sequently, when comparing observational data that span 
a wide range of scales (Alnfcjj 12 in this case) it makes 



sense to account for the variation of n{k) with scale in 
addition to the so-called "tilt" of the spectrum. We ac- 
count for the running of the spectral index by calculating 
it in specific inflationary models and show that the run- 
ning can have an important effect on structure on galac- 
tic scales. We also estimate the effect of a "hot dark 
matter" component in the form of massive neutrinos on 
the central densities of dark matter halos. We make our 
estimates for halo densities using a semi-analytic model 
normalized against N-body simulations. The model relies 
on our, now well-founded, understanding that halo cen- 
tral densities are determined by their collapse histories 

Although we examine models with varying amounts 
of small-scale power, we are not free to alter the spec- 
trum by an arbitrary amount. The amplitude of power 
on small (~ 8 h^^ Mpc ) scales is often quantified in 
terms of ag, which is the rms overdensity smoothed with 
a top- hat filter of radius 8 h~^ Mpc . Observationally, 
this quantity can be determined in a number of ways 
including techniques that rely on the abundance of rich 
x-ray clusters, the cosmic shear from weak gravitational 
lensing, and galactic peculiar velocity flows. However, 
these estimates do not converge on a definitive value 
(even when the same method is used by different authors) 
and many recent estimates seem to advocate surprisingly 
low values of as (2^, |28). Roughly speaking, recent es- 
timates yield values that span the range 0.55^ 1-2 
for 0.3< flM^ 0.5. In the following, we will only consider 
models with erg > 0.55 because models with a smaller erg 
do not have a good chance of being able to match obser- 
vations. We acknowledge that even this limit is pushing 
the observational bounds (although it is consistent with 
Ref. |27|) but we feel that it is best to explore all pos- 
sibilities for the sake of alleviating the tension between 
theory and observation on subgalactic scales. For refer- 
ence, a scale invariant spectrum (i.e., with spectral index 
n ~ I) that is normalized to the Cosmic Background Ex- 
plorer (COBE) measurements of the large-scale cos- 
mic microwave background (CMB) anisotropy via the fit- 
ting forms of Bunn, Liddle, and White has cg ~ 0.95 
assuming that the gravitational wave contribution is neg- 
liglible. 

The remainder of this paper is organized as follows. In 
Sec. II we introduce several models of inflation and cal- 
culate the power spectra predicted by each model. We 
give a short description of the effects of massive neutrinos 
on the evolved linear power spectrum in Sec. HI. We dis- 
cuss the properties of dark matter halos and describe our 
semi-analytic model for estimating halo central densities 
in Sec. In Sec. ^ we present our results and compare 
them with the observed central densities of dwarf and low 
surface brightness (LSB) galaxies. Lastly, we summarize 
our conclusions and indicate directions for future work in 



Sec. VI. 



Throughout this work we will assume flu = 0.3, f^A 
0.7, QbIi^ = 0.020, and h = 0.72. 
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II. INFLATIONARY POWER SPECTRA 

It is widely believed that the primordial density pertur- 
bations that led to the growth of structure in the Universe 
were produced during inflation: quantum fluctuations in 
the inflaton field were frozen in as the rapid cosmological 
expansion stretched these fluctuations to length scales 
larger than the horizon. The power spectrum of primor- 
dial perturbations can be calculated via the slow roll ap- 
proximation (for a review, see Ref. [Sl] and references 
therein). The standard calculation, to lowest order in 
slow roll, yields expressions for the spectrum of density 
perturbations at horizon crossing. 



S^k) 



1 V^(b) 



(1) 



and the effective spectral index of the primordial spec- 
trum, 



n(fc) ~ l + 2?7-6e 



(2) 



in terms of the inflaton potential V{(j)) {V'{(p) = dV/d(p, 
V" = (PV/d4P') and the slow roll parameters e = 
mlJ2{V'/Vf and 77 = mli{V" /V). The reduced Planck 
mass is defined in terms of Newton's constant as nipi = 
I/a/SttGn — 2.4 X 10^* GeV and, as usual, these expres- 
sions are to be evaluated at horizon crossing (i.e., when 
k = aH). 

In the limit of exact de Sitter space during infiation, 
the predicted primordial power spectrum would approach 
exact scale invariance; however, any model in which a 
scalar field is slowly rolling towards a minimum of its po- 
tential will predict some deviation from n = 1. Of course, 
in the context of some inflationary models, the deviation 
is quite small, and scale invariance is a reasonable ap- 
proximation. One frequently cited example of this sort is 
power-law inflation, for which there is an exact solution 
. The reason why approximate scale invariance is ex- 
pected in this model has to do with estimates of the grav- 
itational wave contribution to the CMB quadrupole. In 
addition to scalar density fluctuations, inflation also pro- 
duces tensor (gravity wave) fluctuations. In the power- 
law inflation case, the ratio of the tensor to scalar con- 
tribution to the CMB quadrupole, r = Ctensor / ^scalar ^ 

increases with the tilt as r ~ 6.9(1 — n). A similar result 
also applies to chaotic inflation models [e.g., Ref. |40 ) 
because (p ~ rripi in these models. Recent CMB mea- 
surements indicate that the tensor contribution is small 
(r^ 0.2 ^) so power-law inflation requires 0.97. 

However, this case does not exemplify general infla- 
tionary predictions. Tensor perturbations can be negli- 
gible even if the tilt is not. This is because the tensor 
wave amplitude depends only on the energy scale of in- 
flation. The gravity wave contribution is negligible if 
the inflaton fleld remains far below the Planck scale, as 
would be expected in well-motivated models such as the 
running-mass case discussed below. Models meeting this 



requirement can naturally produce modest tilts and spec- 
tral index running. Moreover, there are reasonably well- 
motivated cases that can yield dramatic departures from 
scale invariance. 

In the balance of this section, we briefly outline the 
predictions of several models of inflation that lead to de- 
viations from the standard n = 1, scale invariant pri- 
mordial spectrum and present the z = linear power 
spectra in each case. Included in this set of models is 
a more extreme example that exhibits so-called "broken 
scale invariance" and for which the slow roll approxima- 
tion cannot be used. In all other cases, we calculate the 
primordial power spectrum to second order in slow roll 
using the method of Stewart and Gong ]35[ | which is suf- 
flciently accurate for our purposes ||36| . In this way, we 
take into account both the tilt of the power spectrum and 
the running of the spectral index. 

To derive the low-redshift power spectra, we use the 
fltting form for the transfer function given by Eisenstein 
and Hu and the exact relation for the linear growth 
factor in flat cosmologies with a cosmological constant 
given by Bildhauer, Buchert, and Kasai [Q. In all cases 
we normalize the power spectrum to COBE using the flt- 
ting functions of Bunn, Liddle and White Q . We con- 
sider several models in which the effective spectral index 
varies significantly with scale. In these cases, wc follow 
the prescription of Ref. poj and evaluate the normaliza- 
tion at the scale fc* = IHq ~ 0.0023 h Mpc~^, which is 
approximately the pivot scale of the COBE data, using 
the effective spectral index at that scale, n{k^). 



A. Inverted power law potentials 

We begin with the illustrative example of the inverted 
power law (IPL) potential (or "small field polynomial" 
in the language of Ref. |Q ) which has the basic charac- 
teristics of "new inflation" p9|]. The general form is 



(3) 



with p > 2. This potential implies that the effective 
spectral index of the primordial power spectrum on the 
scale fc* is given by 



n(fc*) 



1 




(4) 



where is the number of e-folds of inflation that occur 
between the epoch when the scale fc» leaves the horizon 
and the end of inflation. We can obtain a rough esti- 
mate of N^, in terms of the energy density at reheating, 
Pkh, the value of the inflaton potential when fc* leaves 
the horizon, , and the value of the inflaton potential at 
the end of inflation, Vp, by assuming instantaneous tran- 
sitions between vacuum domination and matter domina- 
tion at the end of inflation and matter domination and 
radiation domination at reheating. This gives 
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If the details of the end of inflation and the process of re- 
heating were known, would be known precisely; how- 
ever, these details are not known. In order to obtain 
definite predictions, we take = 50 which is a fairly 
standard working hypothesis. Using this in Eq. (U), we 
see that with p = 4 (we refer to this model as IPL4) 
this model predicts mild deviation from scale invariance, 
namely ~ 0.94. Accordingly, the spectral index is 

mildly scale-dependent, \dn{k)/d\n{k)\ ~ 0.001. Figure 
1 depicts a typical power spectrum at z = with the 
choice p — A. Rather than P(fc) or A^(fc), we plot the 
rms overdensity on a given mass scale cr(M), because this 
is the relevant quantity for our subsequent calculations. 
The COBE normalization amounts to choosing a suitable 
combination of Vq and c and the effective spectral index 
is insensitive to this choice. Normalized to COBE, this 
model predicts a perfectly acceptable value of as — 0.83.^ 
Before proceeding, we mention that the particle 
physics motivation for this type of potential may be 
somewhat dubious. In particular, if p = 4, COBE nor- 
malization requires the dimensionless 0* coupling con- 
stant to be of order 10"-'^'* (the fine-tuning problem may 
be obviated by considering the coupling to be a parame- 
ter of an effective field theory rather than a fundamental 
parameter) Nevertheless, we consider this to be a 

good illustrative example because it is simple and has the 
general behavior \n{k) — 1| ~ 0{1)/N{k) that is exhib- 
ited by a wide variety of models including many specific 
incarnations of new inflation |39[| , hybrid inflation and 
mutated hybrid inflation as well as models with a 
variable gravitational constant . In addition, this po- 
tential mimics the potential encountered in a particular 
variation of mutated hybrid inflation known as "smooth" 
hybrid inflation [||. 

B. Running- mass model 

Stewart has proposed a model in which the need to 
fine-tune the inflaton mass in order for inflation to occur 
in the context of supergravity is eliminated by a "flatten- 
ing" of the effective inflaton potential due to loop correc- 
tions . This provides a natural mechanism for gener- 
ating a potential that gives rise to inflation. Interestingly, 
the resulting effective potential can lead to a spectral in- 
dex considerably different from n = 1 and with a signifi- 
cant scale dependence. In this model it is assumed that 
in the sector of the infiaton field, supersymmetry is bro- 
ken explicitly during infiation and the scalar fields have 



^ The length scale 8 h ^ Mpc corresponds to a mass scale of M ~ 
1.8 X 10^* h~^ Mq under the assumption that = 0.3. 



soft supersymmetry-breaking mass terms as would gen- 
erally be the case. Through couplings to fields with soft 
supersymmetry-breaking masses, the scalar field masses 
may get important renormalization corrections. The one- 
loop correction to the inflaton potential then gives an 
effective potential with a running inflaton mass, 

V{4>)^V^ + Kn'{cj,)<f' + ..., (6) 

where the ellipsis represents nonrenormalizable terms 
that become important at the Planck scale. The value 
of Vo is tied to the scale of the supersymmetry breaking 
during inflation, Afs V^^'^. 

This type of model has been discussed extensively by 
Covi, Lyth and Roszkowski and Covi and Lyth 
who derived cosmic microwave background constraints 
on such models. Inflation occurs in the vicinity of an 
extremum of the potential and the established blueprint 
for analyzing inflation in the context of this model is to 
assume that ■m'^{(p) can be approximated by a function 
that is linear in ln((/!)) while cosmological scales are leaving 
the horizon Q In the notation of Refs. 0, |6|, 

the effective potential is then written as in Eq. (ra) with 
m'^{(f>) ~ — (Vb/mp;)(c/2 — cln((/)/^o)) and the spectral 
index in this case is given by 

n(fc,) ~ 1-f 2cre-'=^* -2c (7) 

where 0o is chosen such that V'{(j)o) — and c and a 
are parameters that may be either positive or negative 
and we generally expect that |c|^ |(t|< 1 With 
c > and a < 0, we have the particularly interesting 
case that n < 1 and decreases with increasing k. With 
c < and cr > 0, n > 1 and decreases with scale. We 
use such a model in several places to illustrate the pre- 
dictions of a model with a primordial power spectrum 
that has n > 1. For reasonable parameter choices there 
may be a significant tilt and the scale dependence may 
be as strong as |(in(fc)/c?ln(A:)| ~ 0.005 on cosmologi- 
cally interesting scales. The significant running of the 
tilt is not surprising. As we stated earlier, in the con- 
text of slow roll inflation, models with significant tilts 
typically exhibit strong variation of n with scale. The 
COBE normalization is related to the parameters Vq and 
r = |c| hi{mpi / (j)o) . In each case we make physically rea- 
sonable choices of these parameters to enforce the COBE 
normalization (see Ref. ) . The shape of the spectrum 
is relatively insensitive to these choices. We choose this 
model as a simple example of an inflationary model that 
may naturally predict significant deviations from n = 1 
and running of the spectral index. 

In Fig. 1 we show present-day, linear power spectra 
for two particular choices of parameters. The model with 
n > 1 (RM n > 1) has parameter choices a — 0.05 and 
c = —0.001, resulting in n(fc*) ~ 1.1 and as ~ 1.21, 
and is shown largely for illustration (note that in this 
case a hybrid mechanism is necessary to end infiation). 
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The more interesting case with n < 1 (RM n < 1) has 
parameter choices a — —0.31 and c ~ 0.04. For this 
model, n(fc*) ~ 0.84 and |dn/(iln(fc)| ~ 0.004 which is 
consistent with constraints on tih from various analyses 
of CMB, large-scale structure and Lya forest data ||3^, 
|33| . For this model as — 0.65 which is on the low side of 
our acceptable range. 

C. Spectra with broken scale invariance 

In contrast with the above models, there may be a 
feature at some scale in the inflationary potential that 
causes the power to drop abruptly. This possibility leads 
us to consider models with so-called "broken scale invari- 
ant" (BSI) spectra. In these models, there is a critical 
scale fee J and for kc and k <ti kc the primordial power 
spectrum has an effective power law index rt ~ 1. How- 
ever, on scales near the critical scale the amplitude of the 
initial density perturbations changes abruptly so that the 
power on scales k > kc can be significantly less than that 
on scales k < k^. This type of spectrum may arise in 
models in which more than one field plays an important 
role in inflation while cosmological scales are leaving the 
horizon [Q but placing the scale k^ in an observationally 
interesting range usually introduces a fine-tuning issue. 

As an idealized case of BSI, Starobinsky derived an 
analytic expression for the primordial power spectrum in 
a model where the inflaton potential has a step discon- 
tinuity in its first derivative (c/. Eq. |^)] [|9|. This is a 
useful model to study for two reasons. First, the primor- 
dial power spectrum can be written in a relatively sim- 
ple closed form. Second, and more material, this model 
exhibits the most rapid drop in power possible in a sin- 
gle field model of inflation |^^. Lesgourgues, Polarski 
and Starobinsky investigated using primordial power 
spectra of this type to explain a feature on scales of about 
125 h~^ Mpc in the galaxy power spectrum measured 
by the Automatic Plate-measuring Machine survey js^ 
(such a feature is not present in more recent determina- 
tions of the galaxy power spectrum) while Kamionkowski 
and Liddle |Q explored the effects of such a primordial 
spectrum on the typical abundance of dwarf satellites. 

In this scenario, the power spectrum of density pertur- 
bations Ap(fc), prior to being modified by causal physical 
processes, is given by the following exact relation [ |49| : 
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FIG. 1: The z = Q rms overdensity as a function of mass 
scale predicted by several models of inflation and normalized 
to COBE. The models depicted here are the running mass 
model with a = 0.05 and c = -0.001 (RM n > 1), the n = 
1 scale invariant Harrison-Zel'dovich spectrum (n = 1), the 
broken-scale invariant model with kc = 0.9 and p = 10 (BSI), 
the inverted power law model with index p = 4 (IPL4), and 
the running mass model with a = —0.31 and c = 0.04 (RM 
n < 1). 

normalization is according to the COBE data. Inspired 
by the work of Kamionkowski and Liddle |^ , we choose 
p = 10 and kc = 0.9 h Mpc^^ in order to suppress power 
on mass scales M"^ 10^" h^^ Mq and thereby alleviate 
the dwarf satellite problem (note that we have chosen a 
different kc than Kamionkowski and Liddle |^6| , partly 
because we assume a different cosmological model) . The 
dotted line in Fig. |^ shows the 2 = 0, linear power spec- 
trum predicted by this model. As there is a rise in power 
prior to the cutoff at M ~ 10^° h"! Mq (see Fig. |), 
we find that erg — 0.97 which is slightly larger than the 
value in a scale invariant model. 
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(8) 



where y = k/kc, /± = 1 ± l/y^, and p is the ratio of the 
amplitude of fluctuations on scales k < kc to the ampli- 
tude of fluctuations on scales k > kc- In this model, we 
are free to choose the amplitude of primordial tensor per- 
turbations and we assume that they are negligible. The 



III. MASSIVE NEUTRINOS 

A preponderance of evidence from solar and at- 
mospheric neutrino oscillation experiments like Super- 
Kamiokande [Q, the Sudbury Neutrino Observatory 
p3| , the Russian- American Gallium Experiment , the 
Gallium Neutrino Observatory |Q , the Gallium Experi- 
ment p6| , and the Soudan Experiment seems to im- 
ply that neutrinos are indeed massive. Yet these exper- 
iments cannot determine the absolute magnitude of the 
neutrino masses and it may be that the masses are large 
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enough to have significant cosmological implications. If 
massive neutrinos (or other "hot dark matter" particles) 
make up a non-negligible portion of the dark matter, the 
effect of their free-streaming will be to reduce power rel- 
ative to the standard model on small length scales. This 
situation is commonly referred to as the cold-|-hot dark 
matter scenario [ p8[ . 

It is easy to estimate the scale at which this effect 
becomes important. Massive neutrinos will move at 
a speed over order c until they become nonrelativistic 
when TOy ~ 3Ti^ which occurs at a redshift of znr — 
2 X 10'^(my/eV) and we expect power to be suppressed 
on scales smaller than the horizon scale at redshift znr. 
As such, a rough estimate is that power will be damped 
on all scales kps where 



fcps - o.omli^^l^h Mpc-1. 



(9) 



This corresponds to suppression of power on mass scales 
M< Mfs =^ 3 X 10i8f2-3/2(m^/eV)-3/2 h"! Mq . The 
contribution of N^, massive, light (m^ <C 1 MeV) neutri- 
nos to the mean matter density, relative to the critical 
density, is fli, ~ N,^{m^ /eV)h~^ /92. On scales k >> 
k-ps, the fractional suppression of power due to massive 
neutrinos approaches a value of ^ (1 4- 8f7y/5lM)^^ [ p9[ . 

This modification to the power spectrum on small 
scales has been studied in detail by many authors 
{e.g., Refs. |37[ ^). In fact, the current best bounds 
on neutrino masses come from demanding consistency of 
the power spectrum on scales probed by COBE and the 
smaller scales probed by clusters and the Lya for- 
est or from the shape of the observed galaxy power 
spectrum on scales 0.01^ k'^ 0.2 |62j. Roughly speak- 
ing, these cosmological bounds dictate that three nearly 
mass-degenerate neutrinos must have TOj/^ 1 — 2 eV |5S|| 
while direct bounds on the electron neutrino mass from 
tritium decay experiments give TOj/^^ 3 eV [ |63[ . 

In what follows, we study the effect of the suppression 
of small-scale power by three neutrinos with effectively 
degenerate masses m,y'^ 1 eV on the central densities of 
dark matter halos. As was pointed out in Ref. pof , ex- 
cessively large neutrino masses lead to unacceptably low 
values of a^. Our strategy is to fix f^M — 1— ^^A = 0.3 and 
n = 1 and to ascertain whether or not a neutrino mass 
that saturates our lower limit of as > 0.55 can alleviate 
the central densities problem associated with the CDM 
paradigm. For this cosmology, the lower limit on as is 
saturated by a neutrino with — 0.65 cV. For compar- 
ison, we also report results for a model with nii, — 0.5 
eV which has as — 0.64. We show in Fig. ||, the present- 
day linear power spectra of the scale invariant reference 
model and the two models with massive neutrinos. No- 
tice the suppression of power on all relevant mass scales. 
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FIG. 2: Power spectra with massive neutrinos compared to 
the standard, scale invariant power spectrum with no massive 
neutrinos. 



IV. THE CENTRAL DENSITIES OF DARK 
MATTER HALOS 

A. Describing dark matter halos 

The absolute size of a virialized dark matter halo 
can be described by the virial mass Mvir, or equiva- 
lently the virial radius i?vir, or the virial velocity V^-^^ = 
GMy-„/ R^i^. The virial radius of a halo is defined as 
the radius within which the mean density is equal to the 
virial overdensity Avir , times the mean matter density of 
the Universe, pM- Thus Mvir and i?vir are related by 



Mvir = Yi?virPMAv 



(10) 



The virial overdensity is set by the spherical top-hat col- 
lapse approximation and for flat cosmologies the value of 
Avir at redshift z can be approximated by ||64| 



Avir(2:) 



IStt^ + 82a; - iQx^ 
x + l ' 



(11) 



where x+l = VIm{z) f^M(l + zf/{nM{l + zf + f^A). 
With Q.M = 0.3 and Q.^ = 0.7, Avir 337 at z = 0. 

Several analytic density profiles for dark matter halos 
have been proposed as good approximations to the results 
of high-resolution N-body simulations and all agree at 
large radii. Moore et al. S found that the density in 
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the central region of the halo varies as p cx r and so 
proposed the density profile 



ir/rMV-^il + r/ru) 



1.5 



(12) 



while Navarro, Frenk, and White (hereafter NFW) 
found that density in the central regions of halos varies 
as p cx r^"'^ and have therefore proposed the profile 



Pnfw(0 



(r/r,)(l +r/r,)2 



(13) 



In most of what follows, we will assume an NFW profile. 
For the NFW profile the two parameters are related by 
Ps ~ pNFw(0.466rs). If N-body simulations do predict 
profiles somewhere between the p cx r^^ and p oc r~^'^, 
the behaviors of the NFW and Moore profiles respec- 
tively, then the NFW profile is a conservative choice in 
the sense that Moore profiles predict more centrally con- 
centrated halos 0. NFW profiles can more easily match 
the data and require less drastic modifications of the 
standard paradigm in order to do so. 

A useful criterion for assessing the relative central con- 
centration of a halo is the concentration parameter 



R 



vir 



Although this quantity is defined explicitly in terms of a 
parameter of the NFW profile, viz. Ts, there is no signif- 
icant loss of generality because we can identify Vg with 
the point at which dlnp/dlnr = —2 and thereby relate 
these results to the results obtained using another den- 
sity profile Q . Restating the NFW profile as a velocity 
curve gives 



V^(r) = V^- — ^ 

* C \ J * vir r / 

/(Cvir 



f{x) 



(15) 



with X = r/vg and f{y) = ln(l + y) — 2//(l + J/)- The 
maximum velocity is given by 



2 

max 



0.216K', 



/(Cv 



(16) 



and occurs at a radius 



2.16rs 



A more directly observable measure of halo central 
densities has been proposed by Alam, Bullock, and Wein- 
berg jj]. This quantity is known as the "central density 
parameter" and is defined as 



half of Vrr, 



In practical units, A 



V/2 



5 X 



10^(14iax/100 kms~^)2(rv/2/h~^ kpc)'^. Assuming an 
NFW profile, rv/2 — O.lSrg and the central density pa- 
rameter is given by 



A^™ 340cL//(Cvir). 



(18) 



There are distinct advantages to using Av/2 to charac- 
terize the central densities of dark halos. For one, Av/2 
is more robustly determined, both observationally and 
in numerical data from N-body simulations, than is the 
inner slope of a density or velocity profile yet it probes 
scales small enough to betray the conflict between the- 
ory and observation. Moreover, Av/2 is defined without 
reference to any particular density or velocity profile. 



B. Predicting halo central densities 

To find the central densities of dark matter halos pre- 
dicted by the aforementioned inflationary models, we 
make use of the semi-analytic model of Bullock et al. |2C[| 
who were stimulated by the previous work of NFW [p5|. 
This model has been calibrated to the results of high- 
resolution N-body simulations and although the model 
was developed in the context of scale invariant CDM 

(14) 

power spectra, ^ it was shown to work well in predict- 
ing the results of an LCDM simulation with n — 0.9, 
as discussed in Ref. [Q. This model represents an im- 
provement over the previous NFW model because it re- 
produces the relationship between Cvir and Mvir observed 
in N-body simulations as a function of redshift whereas 
the NFW model fails at z > 0. It is important to realize 
however, that the following treatment is simplified and 
untested over the full range of power spectra we apply it 
to. Specifically, this model has not been tested against 
simulations with running spectral indices nor has it been 
tested against simulations with a significant hot dark 
matter component. Note, however, that halos formed 
within hot -|- cold dark matter simulations do seem to 
follow an NFW profile down to ~ 2% of the halo virial 
radius |69| . Ideally, our results will motivate future work 
using N-body simulations in order to test our preliminary 
conclusions. 

Briefly, the model of Bullock et al. ^ (and NFW) 
embodies the fact that we expect the central densities 
of dark matter halos to reflect the mean density of the 
Universe at a time when the central region of the halo 
was accreting matter at a high rate ||l^. Therefore, we 
expect halos with central regions that collapsed earlier 



Pcrit 2 y iiorv/2 y 

or the mean dark matter density within the ra- 
dius rv/2 where the galaxy rotation curve goes to 



^ The Bullock et al. model was shown to work well in predicting 
the redshift and mass dependence of halo profiles for a standard 
LCDM model, and also reproduces the z = results presented 
by NFW for standard CDM, open CDM, LCDM, and several 
power-law models |20j. 
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to be denser than their late-forming counterparts. Ac- 
cordingly, our first step is to assign an epoch of collapse 
to a halo via the prescription that, at the collapse epoch 
Zc, the typical collapsing mass, M^[zc), is equal to some 
fixed fraction of the halo's virial mass. Explicitly, we 
define 



(19) 



M*(zc) is the mass scale at which the rms density fluc- 
tuation is equal to the equivalent linear overdensity at 
collapse, (5c — 1.69. If cr(M, z) is the rms overdensity on 
mass scale M at redshift z [we use a{M) with no red- 
shift argument to denote cr(M, z = 0) as usual], then 
the collapse criterion can be written as a{M^, Zc) = 1.69. 
Notice that this definition of the epoch of collapse differs 
from that of NEW who defined the collapse epoch using 
the extended Press-Schechter formalism [Q. This is the 
key difference that gives the improved model the ability 
to trace the redshift dependence of the Mvir — Cvir rela- 
tionship, is a free parameter and Bullock et al. found 
that the model is in good agreement with the results of 
N-body simulations if F = 0.01 [|o). The smaU value 
of the parameter F is not surprising. The densities that 
characterize the very central regions of halos are deter- 
mined by the power on scales much smaller than the size 
of the halo, scales that broke away from the expansion at 
a much earlier time than the mass scale Afvir. 

It is already evident that the central densities of dark 
matter halos will be very sensitive to a{M) on small 
scales and hence, to the slope of the primordial power 
spectrum or the presence of hot dark matter. At early 
times (but during matter domination), CIm ~ 1, and 
a{M,z) oc (1 -t- z)^^. Thus the epoch of collapse varies 
approximately as (1 -I- Zc) oc a{FMvir)- If the central 
densities do, in fact, reflect the mean density of the Uni- 
verse at the epoch of collapse then, roughly speaking, we 
expect Av/2 oc (T'^(FMvir) so that a change in power by 
a factor of 2 will lead to a change in central density by a 
factor - 8 [cf. Eq. (|l|) and Eq. (||) below]. 

The second step is to relate the mass density of the Uni- 
verse at Zc to a characteristic halo density today. Bullock 
et al. chose to use ps defined by 



An 



(20) 



For an NEW profile, ps — 3psf{cvii)- Introducing the free 
parameter K, we associate ps with the universal density 
at collapse via 



Agreement with N-body simulations fixes K = 4.0. Bul- 
lock et al. and Wechsler et al. claim that the la scatter in 
the Cvir — Afvir relation is roughly A log(cvir) — 0.14 while 
Jing has argued for a somewhat smaller scatter given by 
Alog(cvir) ^ 0.08 at la [|7|. 

With this model in place, we can use the linear power 
spectra of the previous two sections to predict Cvir and, 
more practically. Ay /2 and compare them with observa- 
tion. 



V. RESULTS 

In this section we compare the predictions of the pre- 
vious models with data on the rotation curves of dwarf 
galaxies and low surface brightness (LSB) galaxies. We 
concentrate our discussion on galaxies with both HI and 
Ha data or HI data that has been corrected for the effects 
of beam smearing. The data we use are taken from the re- 
cent works of Swaters Q (mass modeling of these galax- 
ies was performed by van den Bosch and Swaters [[ll| ) de 
Blok, McGaugh and Rubin and de Blok and Bosma 
jlOf who combined previous HI measurements |68[ |7^ 
with high-resolution Ho; measurements. We use these 
data to derive observational estimates of Av/2 for com- 
parison with the theoretical predictions. Eor the data of 
de Blok and Bosma [|l^, we use their best fitting model 
for the dark matter distribution of each galaxy in the ab- 
sence of baryons to derive estimates for Av/2. For the 
data of Swaters j6^ and de Blok, McGaugh, and Rubin 
Q, we fit the raw data to the velocity profile proposed 
by Kravtsov et al. ]|69[, 



Vc{r) = V^- 



[l + (r/rfc)"]('r+/3)/"' 



(23) 



and use the best fitting models to estimate Av/2. The 
profile in Eq. ( p3| ) has the practical advantage that it pa- 
rameterizes the sharpness of the transition between the 
two power laws at large and small radii. Hence the fitted 
value of the effective power law index at small radii is to 
some degree decoupled from the details of the rotation 
curve at r^. This added versatility makes it a very 
useful and accurate formula for describing observed ro- 
tation curves at small radii. Our estimates of Av/2 are 
robust in that for most galaxies in the above samples, 
the inferred values of Av/2 change by less than 40% if 
we instead fit the data with NEW, pseudoisothermal or 
Burkert |^ density profiles.^ The robustness of the cen- 
tral density parameter is yet another advantage of using 
Ay /2 as a diagnostic of the central densities of dark mat- 
ter halos. 



Ps = i^VcritAvir(2 = Q)^m(1 + Zcf ■ (21) 

Solving Eqs. ( |20| ) and (^ for Cvir gives 

Cvir(Mvir) = K{1 + Zc(Mvir)). (22) 



It is interesting that there is no systematic difference in the de- 
rived Av/2 from one profile to the next. However, Moore profiles 
tend to fit the data more poorly, and give larger variation in the 
implied Av/2. This is similar to the result found by van den 
Bosch and Swaters 
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Any comparison of the predictions of N-body simula- 
tions or semi-analytic calculations that model the behav- 
ior of CDM with data rests on some assumptions about 
the physics of baryon infall. In making this comparison, 
we believe that our methods are conservative in the sense 
that we likely overestimate Av/2 based on the observa- 
tional data in order to give the data every opportunity 
to match theoretical predictions (including the scale in- 
variant "standard model" ) . First, we restrict ourselves to 
dwarf and LSB galaxies which are generally believed to 
be dominated by their dark matter components [7l| . 
In so doing, we hope that any effects of baryonic infall are 
mitigated but recognize the fact that we may be intro- 
ducing a heretofore unappreciated selection effect. Sec- 
ond, we calculate Av/2 based on the full rotation curve 
without mass modeling or baryon subtraction. We there- 
fore overestimate the central density of the primordial 
dark matter halo because the cooling and contraction of 
the baryons likely lead to contraction of the dark matter 
component as well (see Ref. Q, although also see Ref. 
p7t). Third, the measured rotation curves of about 20% 
of the galaxies in the sample may not extend out to large 
enough radii for an accurate determination of Knax and 
consequently Vmax may be significantly underestimated 
for several galaxies. In these cases, we simply take the 
last point in the rotation curve as an estimate of Knax- 
By examining Eq. (17) it is easy to see that if Vc{r) oc r'^ 
with 7 < 1 at small radii, an underestimation of Knax 
by a factor /v leads to an over estimation of Av/2 by a 

factor of (clearly, for 7=1, corresponding to a 

constant density core, the error cancels exactly). In other 
words, the error introduced has the net effect of bringing 
theory and observation closer together. 

In Fig. ||, we show the theoretical predictions for the 
concentration parameter Cvir in the context of our infla- 
tionary models. Figure ^ shows the predictions in scenar- 
ios with massive neutrinos. Notice the wide swath of the 
Cvir — Afvir plane that is carved out by the various models 
and, in particular, that Cvir can be reduced by a factor 
of two or more by adopting primordial power spectra 
predicted by reasonable models of inflation or by adding 
neutrino masses that are not ruled out by observation or 
experiment. Dark matter halos may be significantly less 
concentrated than standard LCDM plus scale invariance 
predicts. 

Unfortunately, the Cvir — Afvir relation is not directly 
observable and, what is more, it is defined in terms of 
a particular density profile. To connect theory with ob- 
servations, we compare the quantity Av/2, as a measure 
of inner halo concentration, to Knax as a measure of the 
absolute size of the halo. For an NFW proflle, Vmax is 
related to M^-„ through Eqs. ^ and 

The results of this comparison are shown in Figs, ^and 
^. First, consider the predictions of the various models 
of inflation. Although the agreement or disagreement of 
a particular model with the data is hard to quantify, it is 
not surprising that the running mass model with n > 1 
is effectively ruled out by the data. More interestingly. 
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FIG. 3: The median Cvir — Mvir relation predicted by several 
different primordial power spectra. The predictions corre- 
sponding to the diflterent primordial power spectra are labeled 
in the same fashion as in Fig. ^. Bullock et al. ^] estimate 
the la scatter to be Alog(cvir) — 0.14 while Jing argues for a 
smaller scatter of Alog(cvir) — 0.08 These estimates for 
the la scatter are illustrated by the error bars in the upper 
right corner. 



we find that in agreement with previous studies Q, the 
n = 1 scale invariant spectrum also has difficulty repro- 
ducing the observed galactic central densities. This is a 
restatement of the problem: if we are not preferentially 
selecting low density galaxies by restricting attention to 
low surface brightness and dwarf galaxies, then some ad- 
ditional physics is needed to reconcile the standard model 
of CDM plus scale invariant primordial spectrum with 
the observed central densities of dark matter- dominated 
galaxies. IPL4 does a somewhat better job of matching 
the data but the moderate tilt and spectral index running 
in this model are likely not sufficient to bring theory and 
observation together. For BSI, the agreement is much 
better but note that it is difficult to lower the theoretical 
Av/2 values further by adjusting the parameters of the 
model. Increasing p, the ratio of power on scales k « kc 
to power on scales k » kc, does not do much to help 
the BSI model come closer to matching the data because 
the fluctuation amplitude cannot drop quickly enough to 
produce a signiflcant decrease in a{M) on scales of in- 
terest. Meanwhile, we cannot decrease kc much further 
without threatening the success of the standard model 
on large scales. 

Notice that the running-mass model with n < 1 (RM 
n < 1) is a relatively good match to the median of 
the data in the Knax — Av/2 plane (perhaps even un- 
dershooting the median). It is worth noting that this 
agreement has come without the need to saturate our 
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FIG. 4: The median Cvir — Mvir relation in models with 
massive neutrinos. 



lower bounds on spectral tilt from CMB and large scale 
structure {n « 0.9 ± 0.1, see Refs. or our lower 

limit on as- The central densities of dark matter halos 
are very sensitive to the initial power spectrum and it 
seems as though the predicted central densities of dark 
matter halos in a LCDM cosmology may be reduced to 
acceptable levels by invoking simple and well-motivated 
models of inflation with n < I and/or a running spectral 
index. 

Likewise, in the case of massive neutrinos, we see that 
by saturating our lower bound on erg, we can reduce 
the predicted median value of Av/2 to observationally 
acceptable levels. It seems that three massive neutri- 
nos with 0.5 eV^ 0.65 eV can decrease small-scale 
power enough to provide a pretty good match to the val- 
ues of Ay /2 inferred from rotation curve data. 



VI. CONCLUSIONS AND DISCUSSION 

The central density problem is one of several difhcul- 
ties confronting the standard paradigm of structure for- 
mation. In this paper we explored solutions that do not 
invoke uncertain baryonic physics but preserve the cold 
and coUisionless properties of the dark matter. In Sec. 
we showed that models of inflation that predict mod- 



erate, yet observationally acceptable tilts, 0. 



0.9, 



may provide an acceptable solution to the central den- 
sity problem. These tilts are consistent with the lat- 
est constraints from joint analyses of CMB anisotropy, 
large-scale structure and Lya forest data ^ 
Moreover, these tilts can be produced in well-motivated 
models of inflation. In fact, we worked in the context of 
specific models throughout this paper and in so doing. 
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FIG. 5: Vinax vs. Av/2 predictions compared with data. The 
filled triangles are the data points from de Blok, McGaugh, 
and Rubin 0] , the gray squares are derived from the data of de 
Blok and Bosma lid], and the open pentagons are points de- 
rived from the data of Swaters ||68|. The different theoretical 
predictions are labeled in the same way as Fig. ^. The error 
bars in the upper right corner show the expected la scatter 
in the theoretical predictions. The smaller range corresponds 
to the Jing |^ estimate and the larger range corresponds to 
the estimate of Bullock et al. iGOll. 



we were able to take into account the important effect 
of the running of the spectral index. To illustrate the 
importance of the running, we also considered a "tilted" 
model with no spectral running and n ~ 0.84 (the effec- 
tive tilt of the RM n <1 model on the COBE scale) and 
found that this model predicts central densities that are 
more than 40% larger than the those predicted by the 
RM n < 1 model over the range 30 kms~^ < V^ax < 200 
kms^^. 

Given that precise measurements of the tilt of the 
power spectrum and the running of the spectral index us- 
ing the data from the Sloan Digital Sky Survey ||7^ (URL 
http: / /www.sdss.org/ ) and the Microwave Anisot ropy 



Probe (URL |http://map. gsfc.nasa.gov/index. html ) are 
on the horizon ||73|] , it may be useful to adopt an empirical 
stance and consider the maximum tilt and running that 
are acceptable with respect to galactic central densities 
without linking the tilt and running through a particular 
inflationary model. As it is difficult to quantify the agree- 
ment or disagreement of a particular parameter choice 
with the data and because the current data certainly do 
not constrain the slope of the relationship between Knax 
and Av/2; we adopt the somewhat arbitrary, but sensi- 
ble, criterion that a model predicts unacceptably diffuse 
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FIG. 6: Vmax vs. Av/2 predictions in models with massive 
neutrinos compared with data. The data points and the error 
bars in the upper right corner are the same as in Fig. H. 



galaxies if Av/2 < 10^ at Vmax = 100 kms"-^. Using the 
this criterion, we find that a lower limit on n(fc,) allowed 
as a function of dn{k^,) / d\nk can be approximated as 



n{K) + Q.7Mn{K)/d\nk> 0.77 



(24) 



These maximally tilted models have cts > 0.55 for 
n{k^)^ 0.75. If we adopt the criterion that a "good" 



fit to the data has A 



V/2 



3 X 10^ at 



100 kms' 



then a good fit to the data is given approximately by 



n(fc,) + 6.97dn(/c*)/dlnA: ~ 0.87. 



(25) 



We also showed that massive neutrinos with 0.5 eV 
^ TOj^^ 0.65 eV may provide an alternative solution to 
the central density problem; however, we consider this 
solution rather less attractive. In order for neutrinos to 
solve the central density problem, it is necessary to nearly 
saturate our adopted lower limit on erg because, relative 
to the standard scale invariant model, the power spec- 
trum is damped by a factor ~ (1 + / on scales 
smaller than ^ 10^® Mq (corresponding, roughly speak- 
ing, to fc ^ fcps) whereas in the inflationary models, 
power falls off continuously with increasing wave num- 
ber. The range of neutrino masses allowed by the above 
criterion that the dark matter halos not be too diffuse 
is m^"^ 0.9 eV, but as we mentioned earlier a neutrino 
mass greater than '-^ 0.65 eV leads to unacceptably small 
values of (Tg. A neutrino mass of — 0.9 eV implies 
that erg ~ 0.46. 

We did not deal directly with the issue of central slopes. 
The problematic issue here is that cold (and warm) dark 



matter halo densities diverge at small radii whereas galac- 
tic rotation curves seem to be fit better with constant 
density cores |6[ ^ ^, ||, 0. While this is a worri- 
some situation, as we discussed above it is difficult to tell 
the degree to which this is a serious challenge to LCDM. 
First, aZ/ observational errors favor constant density cores 
over cusps. Second, while it has been observed that pseu- 
doisothermal density profiles with constant density cores 
often fit galactic rotation curves better than NFW pro- 
files 1^, ^ |l^ , the conclusion that observations indicate 
halos have cores rather than cusps is a nonsequitur. This 
is because all points on the curve contribute to the fit. 
Rotation curve fits are often largely determined by the 
transition region between the two power laws of the pro- 
file, and may not be faithful representations of the ob- 
served rotation curves at small radii (where there are 
relatively few data points). In addition, van den Bosch 
and Swaters showed that most rotation curves can be ac- 
ceptably fit by divergent density profiles as long as the 
galaxies are much less centrally concentrated than stan- 
dard LCDM predicts Q. To address the inner slope 
issue, it is probably more useful to use a fitting form sim- 
ilar to Eq. (p3|), despite the fact that it is not inspired 
by a theoretical model, because it "decouples" the two 
power laws of the model rotation curve. Our solution to 
the central density problem likely cannot solve the cuspy 
halo issue by itself because central cusps are more a re- 
fiection of the cold and collisionless properties of the dark 
matter than the amount of small-scale power {e.g.^ Ref. 
|]l8| ). Nonetheless, as we have already mentioned, the 
cuspy halo issue is to some degree degenerate with the 
central density problem and it may be that solving the 
latter problem may go a long way toward resolving the 
former. 

A third problem associated with the standard LCDM 
paradigm is the dwarf satellite problem (74|. In essence, 
this problem can be stated in the following way: stan- 
dard LCDM overpredicts the number of satellite halos 
with 10 kms^^^ K„ax~ 50 kms~^ by as much as an order 
of magnitude relative to the number of observed satellite 
galaxies in the local group. As we mentioned earlier, 
Kamionkowski and Liddle investigated solving this 
problem with BSI initial power spectra. It is probable 
that at least part of the solution lies in a feedback mech- 
anism, like reionization suppression |25[] . However, the 
degree of feedback needed will depend crucially on the 
input power spectrum. We examine the subhalo issue 
in the context of inflation in a forthcoming companion 
paper |7^ . Briefly, we find that the discord between the- 
ory and observation can be greatly allayed by considering 
models similar to those studied here and thus, the feed- 
back needed to meet observations can be greatly reduced 
or even eliminated. 

Related to the dwarf satellite problem is the recent re- 
sult of Dalai and Kochanek . The perturbing effect of 
substructure in strong gravitational lenses allowed them 
to constrain the fraction of the host halo mass bound 
up in substructure to be 0.006 < /sat < 0.07 (90% con- 
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fidence). They used this result to limit the tilt of the 
primordial spectrum and put constraints on the neutrino 
mass, and they obtained n > 0.94 and < 0.74 eV 
at 95% confidence. Our results on substructure differ 
from those of Dalai and Kochanek. As we discuss in 
our forthcoming paper |Q, we find that for a host halo 
of the relevant mass, the total mass fraction in subha- 
los is typically larger than the lower limit of Dalai and 
Kochanek (/sat > 0.006) even with significantly tilted 
primordial spectra, 0.8. Thus the tilt of the primor- 
dial power spectrum may not yet be significantly con- 
strained by strong lensing results. However, as we have 
demonstrated here, the long "lever arm" from COBE 
scales to the subgalactic regime offers a potentially use- 
ful avenue for constraining models of inflation. Perhaps 
future lensing studies will provide more significant limits, 
and thus test the intriguing possibility that galaxy rota- 
tion curves are telling us something fundamental about 
the early universe. 
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